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It is well known that there exist many mechanisms that may contribute to neutrinoless double beta 
decay. By exploiting the fact that the associated nuclear matrix elements are target dependent we 
show that, given definite experimental results on a sufficient number of targets, one can determine or 
sufficiently constrain all lepton violating parameters including the mass term. As a specific example 
we show that, assuming the observation of the Oi//?/3-decay in three different nuclei, e.g. ^^Ge, 
^""Mo and ^'"^Te, and just three lepton number violating mechanisms (light and heavy neutrino 
mass mechanisms as well as R-parity breaking SUSY mechanism) being active, there are only four 
different solutions for the lepton violating parameters, provided that they are relatively real. In 
particular, assuming evidence of the Oi^/3/?-decay of ^^Ge, the effective neutrino Majorana mass |m/3/3| 
can be almost uniquely extracted by utilizing other existing constraints (cosmological observations 
and tritium /3-decay experiments). We also point out the possibility that the non-observation of the 
0://3/3-decay for some isotopes could be in agreement with a value of |m^/3| in sub eV region. We 
thus suggest that it is important to have at least two different 0///3/3-decay experiments for a given 
nucleus. 
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I. INTRODUCTION 

After the discoveries of oscillations of atmospheric, so- 
lar and terrestrial neutrinos, one has gained a lot of valu- 
able information regarding the mixing matrix and the 
squared mass differences. The absolute scale of the neu- 
trino mass cannot, however, be determined in such ex- 
periments. Our best hope for settling this important is- 
sue as well as solving a second challenging problem, i.e. 
whether the neutrinos are Majorana or Dirac particles, 
is the observation of neutrinoless double beta decay. 

The total lepton number violating neutrinoless double 
beta decay (Oi^/3/3-decay), 



{A,Z) ^ {A,Z + 2) + 2e- 



(1) 



can take place only if the neutrino is a massive Majo- 
rana particle [l[. The measurement of the Oi^/3/3-decay 
rate could, in principle, determine an absolute scale of 
neutrino mass, solve the neutrino mass hierarchy prob- 
lem and provide information about the Majorana CP- 
violating phases of neutrinos. 

The evidence for a Oj//3/3-decay of ^^Ge has been 
claimed by some authors of the Heidelberg-Moscow col- 
laboration at Laboratori Nazionali del Gran Sasso 
with 



3+0.44 



Ti72(^''Ge) = 2.231-31 



X 10 



25 



y- 
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Such a claim has raised some criticism but none of the 
existing experiments can rule it out 3] . The only certain 
way to confirm or refute this claim is with additional 



sensitive experiments in particular the GERDA ex- 
periment [5|, which plans to start taking data this year. 

There is a general consensus that the Oj^/3/3-decay has 
to be observed at different isotopes. Strong limits on the 
Oj^/3/3-decay half-life have been achieved in NEM03 Q 
and CUORICINO experiments: 



t: 



'^o) > 5.8 X 10"^ y, ri72( Te) > 3.0 x 10"* y. 

(3) 

After neutrino oscillations established that the neutri- 
nos are massive, non-degenerate and strongly admixed, 
naturally most people's attention has focused on the light 
neutrino mass mechanism of the Oi^/3/3-decay. 

It is well known, however, that the Oj^/3/3-decay can be 
triggered by a plethora of other lepton number violating 
(LNV) mechanisms. Among these we should mention the 
exchange of heavy neutrinos, the exchange of SUSY su- 
perpartners with R-parity violating, leptoquarks, right- 
handed W-bosons or Kaluza-Klein excitations, among 
others, which have been discussed in the literature. 

So the day after the Oi//3/3-decay is observed and, hope- 
fully, established in a number of nuclei, the main question 
will be what is the dominant mechanism that triggers the 
decays. 

Possibilities to distinguish at least some of the possible 
mechanisms include the analysis of angular correlations 
between the emitted electrons f§\ , study of the branching 
ratios of Oi^/3/3-decays to ground and excited states 
a comparative study of the Oi^/3/3-decay and neutrinoless 
electron capture with emission of positron (Oi^EC/3"'") [l3| 
and analysis of possible connections with other lepton- 
flavor violating processes (e.g., — 67) [ll|. 
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The main disadvantages of the above approaches are: 
small 0i'/3(3-deca.y rates to excited states, suppressed 
Oi/i5C/3+-decay rates, experimental challenges to observe 
the produced X-rays or Auger electrons and the fact that 
most double /3-decay experiments of the next generation 
are not sensitive to electron tracks. 

In this paper we shall analyze what happens, if sev- 
eral mechanisms are active for the Oi'/J/S-decay. We will 
show that all LNV parameters, including the most inter- 
esting mass term, can be determined provided that Oi//3/3 
data from traditional experiments involving a sufficient 
number of nuclear targets become available. 



II. THE COEXISTENCE OF FEW LNV 
MECHANISMS OF THE Oi/^/3-DECAY 

The subject of interest is a coexistence of the following 
LNV mechanisms of the Oi^/3/3-decay: i) Light neutrino 
mass mechanism, ii) Heavy neutrino mass mechanism. 
Both mechanisms assume only left-hand current weak 
interactions, iii) The trilinear R-parity breaking SUSY 
mechanism generated by gluino exchange. For the sake 
of simplicity we shall assume that the lepton violating 
parameters are relatively real as, e.g., is the the situation 
in the case CP conservation. 

The inverse value of the Oi^/3 /3-decay half-life for a given 
isotope (A, Z) can be written as 



1 



^1/2 



(4) 



Here, 'i]u,N,fip and M^*^ ^ are, respectively, the LNV pa- 
rameters and the nuclear matrix elements (NMEs), in 
the order given above. Each of the NMEs depends, in 
general, quite differently on the nuclear structure of the 
particular isotopes {A, Z), {A, Z+1) and {A, Z+2) under 
study. 

G'^^{Eo, Z) is the known phase-space factor {Eq is the 
energy release), which include fourth power of axial- 
coupling constant gA — 1.25. The G'^'^{Eo, Z) contain 
the inverse square of the nuclear radius R^^, compen- 
sated by the factor R in M"". The assumed value of the 
nuclear radius is i? = roA-^^^ with ro — 1.1 fm. The 
phase-space factors are tabulated in Ref. J^]. 

The lepton number violating mechanisms of interest 
together with corresponding nuclear matrix elements are 
presented briefly below. 



Under the assumption of the mixing of three light massive 
Majorana neutrinos the effective Majorana neutrino mass 
{m/3p) takes the form 



E 



2cCP 



rrii 



(ah rrii > 0) 



(6) 



where Uei is the first row of the neutrino mixing matrix 
and ^ are unknown Majorana CP phases. Wi is the 
light neutrino mass {rrii < 1 eV, i=l, 2, 3). In this case 
only left-handed weak interaction is taken into account. 

The nuclear matrix element M^'^ consists of Fermi, 
Gamow- Teller and tensor parts as 



M. 



Of 



F{u) 



9\ 



M, 



GT(u) 



(7) 



Here, qa is axial- vector coupling constant. The Fermi, 
Gamow- Teller and tensor operators are defined in the 
usual way (see Eq. (|10l) below) with exchange potentials 
as given elsewhere |13| 



B. Heavy Majorana neutrino exchange mechanisms 

We assume that the neutrino mass spectrum include 
heavy Majorana states N with masses Mk much larger 
than the energy scale of the Oz//3/3-decay, Mk ^ 1 GeV. 
These heavy states can mediate this process as the pre- 
vious light neutrino exchange mechanism. The difference 
is that the neutrino propagators in the present case can 
be contracted to points and, therefore, the corresponding 
effective transition operators are local unlike in the light 
neutrino exchange mechanism with long range internu- 
cleon interactions. 

The corresponding LNV parameter is given by 



E 

fc=4 



(8) 



Here, 



is the mass of proton. Uek are elements of 



the neutrino mixing matrix associated with left-handed 
current interactions, ^'j, are CP violating phases. 

Separating the Fermi (F), Gamow- Teller (GT) and the 
tensor (T) contributions we write down 



Ml 



M 



F{N) 



9\ 



M, 



GT{N) 



Mr 



T(N) 



= (o.+iE-.^' 



A. Light Majorana neutrino exchange mechanism 

In the case of light-neutrino mass mechanism of the 
Oi^/3/3-decay we have 



rUe 



(5) 



+ H^^^^\rki)<Jki - H^f\ni)Ski] |0+), (9) 



where 



Ski = H^k ■ iki)im ■ rki) ~ (Tki, aki^ak-ai. (10) 

The radial parts of the exchange potentials can be found 
elsewhere 
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C. R-parity breaking SUSY mechanism 

In the SUSY models with R-parity non-conservation 
one encounters LNV couphngs which may also trigger 
the Oi'PfS decay. Recall, that R-parity is a multiplicative 
quantum number defined by i? = (_i)2S+3B+l (s,B,L 
are spin, baryon and lepton number) . Ordinary particles 
have R = while their superpartners R = —1. The 
LNV couplings emerge in this class of SUSY models from 
the R-parity breaking part of the superpotential 

Wj^^ = K,uUL,El + K,kUQ,Dl + ^i,UH2, (11) 

where L, Q stand for lepton and quark SU{2)l doublet 
left-handed superfields, while E'^, for lepton and down 
quark singlet superfields. This results in a lepton vio- 
lating parameter entering the neutrinoless double beta 
decay, ry^^. 

For simplicity we concentrate below on the trilinear 
A' couplings and write rjp^ = ijy . Under reasonable as- 
sumptions the gluino exchange dominates [3]. We have 



i'2 



6 



(12) 



Here, Gf is the Fermi constant, as — 5|/(47r) is SU(3)c 
gauge coupling constant, rriu^ , mj^ and rrig are masses 
of the u-squark, d-squark and gluino, respectively. 

At the hadron level we assume dominance of the pion- 
exchange mode [13, [l^ . We denote the Oi//3/3-decay nu- 
clear matrix element, M'p' , of Eq. ^ as Ai'} 



ilf with 



CA 



M4 



with c. 



(13) 

= m'^/{mpme) {mA = 850 MeV). The structure 
coefficients of the one-pion a^^ and two-pion mode a^'^ 
are [13, El: a^'' = -0.044 and a^"" = 0.20. The partial 
nuclear matrix elements of the l^p SUSY mechanism for 
the Oi^/3/3-decay process are: 

M^T = (0;i E ^k^i" H^rirki) cj. ■ a„ |0+), 

= (0;i E ^rirki) Ski |0+) (14) 

k^l 

with the radial functions given elsewhere [l3,[il. Under 
these assumptions the obtained nuclear matrix elements 
are given in table U 

In obtaining the nuclear matrix elements we used 
the Self-consistent Renormalized Quasiparticle Random 
Phase Approximation (SRQRPA) [l^ to calculate nu- 
clear matrix elements (NMEs) M^" , and Af^r. The 
SRQRPA takes into account the Pauli exclusion prin- 
ciple and conserves the mean particle number in corre- 
lated ground state. For A=76 and 100 nuclear systems 
the single-particle model space consists of — 5hu! os- 
cillator shells and for A=130 nuclear system we used 



TABLE I. The Oz^/3/3-decay NMEs calculated with the Selfcon- 
sistent Renormalized Quasiparticle Random Phase Approx- 
imation (SRQRPA). The Coupled Cluster Method (CCM) 
short-range correlations calculated with CD-Bonn potential 
are taken into account 17]. g.A ~ 1.25. 
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■'^Ge^'"^Se 


7.98 X 10 


-ib 


5.82 


412. 


596. 


iooMo->i««Ru 


5.73 X 10 


-14 


5.15 


404. 


589. 


i30Te->i»°Xe 


5.54 X 10 


-14 


4.70 


385. 


540. 



— bhuj shells plus Oiii/2 and 0ii3/2 levels both for pro- 
tons and neutrons. In the calculation of the Oi^ 13/3 -decay 
NMEs the two-nucleon short-range correlations derived 
from same potential as residual interactions, namely from 
the CD-Bonn potential [17], are considered. The calcu- 
lated NMEs are given in Table U 



III. CALCULATION AND DISCUSSION 

A. Dominance of a single Oi^/3/3-decay mechanism 

Commonly, it is assumed that a single LNV mechanism 
is responsible for the Oi^/3/3-decay. Let suppose it is the 
light neutrino mass (m^^) or heavy neutrino mass (fyw) 
mechanism. Then, the Oi//3/3-decay half-lives Ti and T2 
for two nuclear systems are related with equation 



\mpf!\ 



1 



\M^\VT~G~ 



ml 



(15) 



Here, Gi is the kinematical factor, while and M^"^ 
are nuclear matrix elements associated with and r/N 
parameters, respectively for the target i. 



B. Two active Oi^/3/3-decay mechanisms 

We will now move into the case of two competing {)v(3f3- 
decay mechanisms representing by the LNV parameters 
mpp and rj, 77 could be 77 at or 77^^. In this case we have 
four different sets of two linear equations: 



±1 



"1/3/3 



±1 



VT2G2 



(16) 

For the absolute value of the LNV parameters we find 
two different solutions, 



± 



(17) 
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FIG. 1. The effective Majorana mass of neutrinos in the case 
of two active mechanisms of the Oi^/3/3-decay, namely hght 
and heavy neutrino exchange mechanisms, as function of pa- 
rameter C (see Eq. (O). T^J^C^Ge) = 2.23 x 10^^ y is 
assumed. Solutions 1 and 2 were obtained for equal and oppo- 
site signs on the left hand side of Eqs. ([16} , respectively. The 
bold point indicates the value of mpp, if the light neutrino 
exchange is the only active mechanism. 
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We note, however, that for 77 = Eqs. P7)) and (IT51) are 
reduced to Eq. P^ . 

By assuming now rj = rj^ the solutions for will 
be analyzed for two different combinations of nuclear sys- 
tems, namely with A = 76 and 100 (case I) and A — 
76 and 130 (case II). An additional assumption is that 
the Oi//3/3-decay half-life of ^^Ge has been measured with 
T^J^C^Ge) given in (g]). In Fig. [T]the two solutions for 

|m^^| are plotted as function of ^, where 



(19) 



Indices 1 and 2 denote the above cases I and II respec- 
tively. The parameter ^ represents the unknown half-life 
of the Oi'^/J-decay of ^°°Mo or ^^OTe. We note that for 
^ = 1 the solution for active only light neutrino mass 
mechanism given by Eq. (jlSp is reproduced and that 
^ = means non-observation of the 0:^/3/3-decay for a 
considered isotope. 

By glancing the Fig. [1] we see that for both combi- 
nations of nuclear systems the two solutions for |?tt.^^| 
exhibit similar behavior. For ^ w 1.1 one of the solution 
can be even equal to zero. The second solution predicts 
> 1 e:V. The current lower limits on the half-life 
of the 0iJl3l3-deca,y of i°°Mo or "°Te restrict the effec- 
tive mass of Majorana neutrinos to intervals (0.0 eV, 
7.2 eV) and (0.0 eV, 3.2 eV), respectively. So, within 
the considered assumptions the claim of evidence of the 
Oj^/3/3-decay of ^^Ge 2] is compatible also with inverted 
(jUi < 50 meV) or normal (rrii w few meV) hierarchy 
of neutrino masses. From Fig. [1] it follows that a small 
improvement of the current half-life limit for ^'^''Te up to 



FIG. 2. The effective Majorana mass of neutrinos in the case 
of three active mechanisms of the Oi//3/3-decay, namely light 
and heavy neutrino exchange mechanisms and R-parity break- 
ing SUSY mechanism with gluino exchange, as function of pa- 
rameter ^Tc. Nuclear systems with A = 76, 100 and 130 are 
considered. T{y'2(^''Ge) is the same as in Fig. [T] ^mo ~ 2.6 
and ^Mo ~ 1.2 are assumed in upper (a) and lower (b) pan- 
els, respectively. For each solution there are given in brackets 
signs in front of term 1/VtG in Eqs. ((20} for ^^Ge, ^°°Mo 
and "«re. 



value 4.1 x 10^^ y (^to — 1-1) would exclude these pos- 
sibilities. Another finding is that the non-observation of 
the Oz//3/3-decay for ^^Mo or ^^Oxe (i.e., C = 0) cannot 
rule out the claim for evidence of the Oiyf3/3-deca,y of ""^Ge. 
This can only happen if, in a more sensitive Ge exper- 
iment like GERDA or Majorana, no Oi^/3/3-decay signal 
will be registered. 



C. Three active Oz^/3/3-decay mechanisms 

In the case of three active Oi'/J/S-decay mechanisms rep- 
resented by the LNV parameters m^^, rjN and rjxi assum- 
ing the measurement of the life time of the Oi^/3/3-decay of 
three isotopes one obtains a set of three linear equations: 



±1 



vTiGi 



+ tjnM^ + r^x'M^' , i = 1, 2, 3 (20) 



The equations (pn|) admit a set of four different so- 
lutions, which are exhibited in Fig. [2l The upper and 
lower panels correspond to ^mo = 2.6 (current limit) and 
^Mo = 1-2, respectively. The allowed ranges of mpp are 
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as follows: i) (1.3,2.2) eV and (5.2,6.1) eV for upper 
panel, ii) (0.0,1.7) eV and (3.2,4.1) eV for lower panel. 
We see that for a given value of ^tc allowed intervals for 
\mfjj3\ depend strongly on the value of ^mo- The upper 
two solutions determining the second interval are already 
excluded by the Mainz and Troitsk tritium experiments 
p^ . In case the claim of evidence of the Oi^/3/3 decay 
of ^^Ge would be ruled out by other experiments, i.e. 
for larger value of T^j2C^^^)^ ^^^Y would decrease and 
might be important. 



IV. CONCLUSIONS 

It has been shown that the extraction of the most 
important neutrino mass contribution, entering neutri- 
noless double beta decay, can be disentangled from the 
other mechanisms, if and when the decay rates in a suf- 
ficient number of nuclear targets become available. In 
the present calculation, to simplify the exposition, we re- 
stricted ourselves in the special case of no right handed 
currents and made the assumption that the LNV parame- 
ters are relatively real. To be more specific, in addition to 
the standard light neutrino mass mechanism of the QvjSP- 
decay, we considered two additional LNV mechanisms, 
namely those involving the exchange of heavy neutrinos 
and R-parity breaking SUSY with gluino exchange. We 
find that this improved analysis leads to completely dif- 
ferent results compared to those of one mechanism at a 
time. It is now possible that larger values of I'ti/j^I can 
be consistent with the data, since the contribution of the 
other mechanisms could be interfering with it destruc- 



tively. 

We specifically discussed the extracted value of the 
effective Majorana neutrino mass , mpp, assuming the 
claim of evidence of the Oi//3/3-decay of ''^Ge [1] as a func- 
tion of half-life data for the two promising nuclei, (^""Mo 
and ^^°Te). We showed that in an analysis including two 
and three nuclear systems there are 2 and 4 different pos- 
sible solutions for |r7i^^|, respectively. One of the solu- 
tions leads to small values of 1^1,3^1, when all mechanisms 
add up coherently. This is compatible also with inverted 
{rrii < 50 meV) or normal (m^ w few meV) hierarchy of 
neutrino masses. Other solutions, however, allow quite 
large values of |?ti^^|, even larger than 1 eV. These can, 
of course, be excluded by cosmology and tritium /3-decay 
experiments. It may not, however, be possible to exclude 
these solutions, if the claim of evidence for ^^Ge would 
be ruled out by future experiments, since, then, the val- 
ues we obtain become smaller than those of the other 
experiments. 

It is thus important that experiments involving as 
many different targets as possible be pursued. Further- 
more, in the presence of interference between the various 
mechanisms, the availability of reliable nuclear matrix 
elements becomes more imperative. 

V. ACKNOWLEDGMENTS 

This work was supported in part by the DFG 
project 436 SLK 17/298, the Transregio Project TR27 
" Neutrinos and Beyond" , by the GraduiertenkoUeg 
GRK683, by the VEGA Grant agency under the con- 
tract No. 1/0249/03 and by a Humboldt Research Award. 
Partial support was also provided by the EU network 
MRTN-CT-2006-035683 (UniverseNet) . 



[1] J. Schechter and J.W.F. Valle, Phys. Rev. D 25, 2951 
(1982); M. Hirsch, H.V. Klapdor-Kleingrothaus, S. Ko- 
valenko, Phys. Lett. B 372, 8 (1996). 

[2] H. V. Klapdor-Kleingrothaus, I. V. Krivosheina, Mod. 
Phys. Lett. A 21, 1547 (2006); H. V. Klapdor- 
Kleingrothaus, A. Dietz, H. L. Harney, I. V. Krivosheina, 
Mod. Phys. Lett. A 16, 2409 (2001). 

[3] A. Faessler, G.L. Fogh, E. Lisi, V. Rodin, F. Simkovic, 
Phys. Rev. D 79, 053001 (2009). 

[4] F. T. Avignone, S. R. Elliott, and J. Engel, Rev. Mod. 
Phys. 80, 481 (2008). 

[5] I. Abt et al. (GERDA Collaboration), 
arXiv:hep-ex/0404039 

[6] R. Arnold et al. (NEMO Collaboration), Phys. Rev. Lett. 
95, 182302 (2005). 

[7] C. Arnaboldi et al. (CUORE Collaboration), Phys. Lett. 
B 584, 260 (2004). 

[8] M. Doi, T. Kotani, H. Nishiura and E. Takasugi, Prog. 
Theor. Phys. 69, 602 (1983); F. Deppisch, Ch. Jackson, 
I. Nasteva, S. Soldner-Rembold, Prog. Part. Nucl. Phys. 
64, 278 (2010). 

[9] S. M. Bilenky, J.A. Grifols, Phys. Lett. B 550, 154 
(2002); F. Simkovic, A. Faessler, Prog. Part. Nucl. Phys. 



48, 201 (2002); F. Deppisch, H. Pas, Phys. Rev. Lett. 
98, 232501 (2007). 

[10] M. Hirsch, K. Muto, T. Oda and H. V. Klapdor- 
Kleingrothaus, Z. Phys. A 347, 151 (1994). 

[11] V. Cirigliano, A. Kurylov, M. J. Ramsey-Musolf and P. 
Vogel, Phys. Rev. Lett. 93, 231802 (2004). 

[12] F. Simkovic, G. Pantis, J.D. Vergados, and A. Faessler, 
Phys. Rev. C 60, 055502 (1999). 

[13] F. Simkovic, A. Faessler, V.A. Rodin, P. Vogel, and J. 
Engel, Phys. Rev. C 77, 045503 (2008). 

[14] A. Wodecki, W.A. Kaminski, F. Simkovic, Phys. Rev. D 
60, 115007 (1999). 

[15] A. Faessler, S. Kovalenko, F. Simkovic, and J. Schwieger, 
Phys. Rev. Lett. 78, 183 (1997); A. Faessler, S. Ko- 
valenko, and F. Simkovic, Phys. Rev. D 58, 115004 
(1998). 

[16] D.S. Delion, J. Dukelsky, and P. Schuck, Phys. Rev. C 
55, 2340 (1997); F. Krmpotic et al., Nucl. Phys. A 637, 
295 (1998). 

[17] F. Simkovic, A. Faessler, H. Miither, V. Rodin, and M. 

Stauf, Phys. Rev. C 79, 055501 (2009). 
[18] E. W. Otten and C. Weinheimer, Rep. Prog. Phys. 71, 

086201 (2008). 



